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Abstract. The screening of a 2p core-hole in Na clusters is investigated using density 
functional theory applied to an extended jellium model with an all-electron atom in the 
center. The study is related to recent experiments at the free electron laser at DESY 
in which photoelectron spectra from mass-selected, core-shell ionized metal clusters 
have been recorded. Relaxed and unrelaxed binding energies as well as Kohn-Sham 
orbital energies are calculated in Perdew-Zunger self-interaction-corrected exchange- 
only local spin-density approximation for valence and 2p core electrons in Na clusters 
up to 58 atoms. The relaxed binding energies follow approximately the metal-sphere 
behavior. The same behavior is seen in the experiment for sufficiently big clusters, 
indicating perfect screening and that the relaxation energy due to screening goes to 
the photoelectron. Instead, calculating the kinetic energy of the photoelectrons using 
unrelaxed binding energies or Kohn-Sham orbital energies yields wrong results for 
core-shell electrons. The screening dynamics are investigated using time-dependent 
density functional theory. It is shown that screening occurs on two time scales, a 
core-shell-dependent inner-atomic and an inter-atomic valence electron time scale. In 
the case of Na 2p ionization the remaining electrons in the 2p shell screen within 
tens of attoseconds while the screening due to cluster valence electrons occurs within 
several hundreds of attoseconds. The screening time-scales may be compared to the 
photon energy and cluster size-dependent escape times of the photoelectron in order 
to estimate whether the photoelectron is capable of picking up the relaxation energy 
or whether the residual system is left in an excited state. 
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1. Introduction 

Core-holes are created when matter is irradiated with energetic photons. The net 
binding energy is obtained from the kinetic energy of the photoelectron upon subtracting 
the photon energy. As the target size increases from the isolated atom via molecules 
and clusters to, ultimately, bulk matter, the shift of the net binding energy contains 
information about the screening capability of the other electrons in the system and 
allows to investigate, e.g., metal-to- nonmetal transitions. In fact, x-ray photoelectron 
spectroscopy (XPS) is well-established and work has been extensively devoted to 
measure core-level binding energy shifts in solids and at surfaces p^. However, the 
systematic study of core-level shifts as a function of the target size became possible only 
recently because of the necessity of having size-selected clusters and tunable, powerful 
x-ray sources such as free-electron lasers available. 

From the theoretical point of view, the calculation of core-level binding energies 
in many-electron systems is non-trivial both conceptually and computationally [H [2]. 
The study of atoms, clusters or bulk, in which suddenly a core-electron is removed by 
single-photon ionization, has a long history, starting probably with Slater's transition 
state theory [3l H] and various kinds of self-consistent field methods (ASCFs) [5l E]. 
ASCF-approaches are based on energy balances between electronic configurations, and 
it is not always clear whether relaxed or unrelaxed configurations should be taken. Here, 
"relaxed" means that the other electrons are allowed to assume the new, energetically 
most favorable configuration — albeit without filling the core-hole. "Unrelaxed" means 
that the other electrons are kept frozen. It is intuitively clear that the unrelaxed 
energy difference between final and initial configuration should be considered if the 
photoelectron escapes too rapidly to notice the relaxation dynamics of the other 
electrons [3 El [9]. 

Methodologically, self-consistent field methods like Hartree-Fock (HF) or density 
functional theory (DFT) (see, e.g., |T0]) are usually employed for the calculation of 
energy differences between final and initial configurations in ASCF approaches. As 
both methods introduce orbitals and orbital energies, the interesting question arises 
how these orbital energies are related to the (relaxed or unrelaxed) binding energies 
pTl El [12]. Clearly, the energy difference depends on initial and final state while the 
orbital energies do not, unless fractional occupation numbers, as in Slater's transition 
state theory, are introduced. 

We use an extended Na-cluster jellium model with an all-electron Na atom in the 
center in order to have core-levels in the system at all. This model may be viewed 
as a hybrid of the pure cluster jellium model [13] and bulk-models in which atoms are 
immersed in homogeneous electron gas O [T3] . We employ DFT to calculate relaxed and 
unrelaxed energy differences between the initial ground state and the final configuration 
with either a valence electron or a 2p-core electron removed. The results as a function 
of the cluster size (up to 58 atoms) are compared with the prediction of the metal 
sphere-model. Time-dependent DFT (TDDFT) (see [15] for a state-of-the-art account) 
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is used to investigate the screening dynamics after the sudden removal of a 2p electron. 
We apply our model in this first study to Na clusters (instead of Pb~-clusters in the 
experiments fi6\ [T7] ) in order to keep the number of inner electrons in the embedded 
atom and thus the numerical effort managable. 

The paper is organized as follows. The basic methodology and theoretical 
background is introduced in section |2l The DFT results on the size-dependence of 
binding energies and the TDDFT results on the screening dynamics are presented in 
section |3l We conclude in section HI 

Atomic units h = = \e\ = ineo = 1 are used, unless noted otherwise. 

2. Basic theory and models 

We use DFT in exchange-only local spin-density approximation (xLSD) with Perdew- 
Zunger (PZ) self-interaction correction (SIC) [H]. The electronic Kohn-Sham (KS) 
equation reads 

eia^^a = (t + V (r) + Vu + Kca " {Vuia + Kci<x}) Via- (l) 

Here, ejo- is the energy of KS orbital ipi„, i = 1,2, . . . N„ with N^- the number of electrons 
with spin-projection cr ='|", J,. T is the kinetic energy operator, V{r) is the spherically 
symmetrical external potential introduced below. Vn and V^ca are the Hartree potential 
and the exchange-correlation (xc) potential in xLSD approximation, respectively (see, 
e.g., [10]). Both are functionals of the (spin) density. The term in the curly bracket is 
the PZ-SIC, i.e., in the KS equation for orbital ipi„ the Hartree-xc potential Vuia + Vxcia 
evaluated with the corresponding spin density Ui^j = |v?jcrp is subtracted from the full 
Hartree-xc potential Vh + Kcco- evaluated with 

n = ^n^ = + n^^, n„ = '^ni^. (2) 

a i=l 

It is known that the PZ-SIC applied to xLSD significantly improves the results as 
compared to uncorrected xLSD [TTl [T9]. The central-field approximation is applied to 
Vr + KcCT — {Vuia + Vxcia} dl]) SO that the single-particle orbital angular momentum 
quantum numbers i and magnetic quantum numbers m remain good quantum numbers. 
Although the PZ-SIC in general improves the uncorrected KS results quantitatively, 
there have been several critical issues discussed in the literature [HI [19], one of them 
being the orbital-dependent effective Hamiltonian in ([T]), leading to non-orthogonal KS 
orbitals. In our calculations we therefore force the KS orbitals to be orthogonal. An 
extended version of the Qprop code [201 121] is used to solve ([T]) and — for the dynamics in 
section [3?2] — its time- dependent version (in which ej^r is replaced by idf and the respective 
adiabatically time-dependent potentials are used [15]). 
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Figure 1. Extended jellium model for a Nasg cluster, (a) Total density (strong 
increase in center because of explicitly treated all-electron Na atom there), (b) KS 
potentials for the various PZ-KS orbitals (gray to black), 1/r dotted; horizontal lines 
indicate KS levels (central atomic ones solid, cluster levels dashed). 



2.1. Cluster jellium model with central all-electron atom 
The external potential V{r) in ([T]) in this work is 



6 — — tor r < K 



for r > R 

The first term —Z/r is the pure Coulomb potential of the central atom for which all 
electrons are taken into account. This atom provides the ionizable core-shell. The other 
A — 1 atoms in the cluster are treated in jellium approximation. The jellium potential 
is that of a homogeneously charged sphere of radius 

R = A^/Vs, (4) 

mimicking the net ionic background of the cluster, as "seen" by the cluster valence 
electrons. A Wigner-Seitz radius for Na of rg = 4 was used. For A = 1 the potential ([3]) 
clearly reduces to the single-atom case. 

Figure [1] shows the results for electron density, the PZ-KS potentials, and the KS 
levels for the neutral closed-shell Nass cluster. Because all electrons in the central atom 
are considered there is a peak in the electron density in the center. What would be the 3s 
level in the isolated Na atom is now the Is jellium cluster level. The higher the jellium 
level the closer the respective KS eigenvalue is to the one in the corresponding pure 
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valence electron jellium model. Note that the PZ-SIC yields the correct 1/r-behavior of 
the KS potentials, unlike pure xLSD. 

2.2. Determination of the core-shell binding energies 

The experimental observable is the kinetic energy ekin of photoelectrons which are 
emitted upon the irradiation of clusters with x-rays of photon energy hw. Energy 
conservation requires 

Eo + nuj = E{ + ekin, ^0 < < (5) 

so that the binding energy can be determined as 

AE = Ef-Eo = huj- ekin, AE > 0. (6) 

Eq < is the initial energy of the target. In what follows we assume that the target 
is initially in its ground state. However, it is far from obvious what Ef is. As long as 
Auger decay occurs on time-scales longer than the photoelectron escape times the final 
target configuration of interest to us is a cluster in which one atom has a core- holeU 
We assume that the photon energy is such that most likely one of the 2p electrons is 
removed. If the escape time is longer than the screening time the photoelectron may 
pick up the relaxation energy 

S = -Eunrlxd — -f'rlxd > 0, (7) 

and the net binding energy is 

AE'rlxd = -E'rlxd — Eq = -Eunrlxd — Eq — T,. (8) 

On the other hand, if the photoelectron escapes before the system is relaxed, the binding 
energy reads 

A-E'unrlxd = -f'unrlxd — Eq = -Erlxd — Eq -\- > AEj-i^^- (9) 

According to Koopmans' theorem the unrelaxed energy difference equals the 
Hartree-Fock orbital energy, — e^^ = -Eun^xd ~ > 0- The fact that in the ultra- 
high photon energy-limit the centroid of the photoelectron spectrum is located at 
Ckin = + hu is known as the Manne-Aberg theorem or the Lundqvist sum rule 
[8]. However, Koopmans' theorem does not hold in DFT. Instead, the "non-Koopmans 
energy" tends to cancel a part of the relaxation energy ([7]) [llj so that KS orbital energies 
|e| = — e are shifted towards AEj-i^d- Whether — e is above or below Ai^rixd depends on 
the xc energy functional used. Introducing fractional occupations < / < 1, Janak's 
theorem establishes a connection between the relaxed binding energy and the relaxed KS 
orbital energies e(/) (see, e.g., [H]). 

Using the extended cluster jellium model introduced in section 12.11 we determined 
the binding energy 

AE,^ = Ei„ - Eq (10) 

X Auger decay will only later lead to secondary, low-energy electrons which, in fact, are seen in 
experiments [1]. 
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of a core-shell electron ia by switching the fractional occupation number of the respective 
KS orbital /jo- from 1 to 0. Without relaxing the final configuration with the core-hole 
we obtain 

^-E'io-.unrlxd = -E'icr.unrlxd ^ Eq, (11) 

and with relaxation 

,unrlxd' 

(12) 

The relaxation energy 

^icr -^jcr, unrlxd -^i(T,rlxd ^-^itr, unrlxd ^-^itr, rlxd; (l*^) 

in general, depends on the orbital from which the electron is removed. Clearly, in the 
case of the removal of a valence electron the relaxation energy is expected to be very 
small because mainly the electrons above the level from which the electron is removed 
contribute to screening. On the other hand, for the screening by the cluster valence 
electrons only the net positive, almost point-like hole of charge |e| counts, and it should 
not matter much from which inner-atomic shell a core-electron is removed. Possible 
differences in the relaxation energy are expected to be due to other core-shell electrons. 
We call this inner-atomic screening, in contrast to inter-atomic screening by cluster 
valence electrons. 

A widely used approach in the DFT community is to use the KS orbital energies as 
a zeroth-order approximation for the binding energy of electrons. However, the initial 
state KS orbital energy may be above or below Ai^jo- ^ixd, depending on the xc potential 
used 111]. Energy differences AE are less sensitive and, even more so, binding energy 
shifts, i.e., the difference of energy differences. 



2.3. Metal sphere model 

The energy required to remove one electron from a metal sphere of radius R which, 
initially, has Zg negative elementary excess charges — |e| reads 

1 — 2Z 

AE^etsph.z. = -Ez^ = W + (14) 

Here, W is the work function. In the bulk limit i? — )■ oo one has Aii^metsph = W, 
independent of Zg. Equation (JH]) can be simply derived by calculating the initial and 
final energies. As the interior of the metal sphere is, by definition, field-free one has 



(with E(r) the electric field) so that 
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Figure 2. Qualitative sketch of the metal sphere model result eq. (fT4|) for the binding 
energy of cluster anions, neutral clusters, and cluster cations vs the inverse cluster 
radius. Starting at — J- (bulk) with the work function W the slope is negative for 
electron removal from anions (dark red, — ), positive for neutrals (blue), and steeper 
positive for cations (green, +). Possible single, isolated atom or ion results are indicated 
by filled circles, which, in general do not lie on the metal sphere lines. The transition 
from the perfect metal behavior to the single, isolated atom or ion result occurs in the 
hatched region. 



The work W needed to remove an electron from the bulk is added "by hand," leading 
to ^\M . Improved metal sphere models yield also W and helped to resolve the "image 
charge paradox" |22j§l . 

Figure [2] summarizes schematically the prediction of the metal sphere model. 
Starting from the bulk value W for the binding energy for — > eq. (031) predicts 
a linear behavior m 1/ R with a negative slope for initially negatively charged clusters, 
a positive slope for initially neutral clusters, and increasingly steep positive slopes for 
more and more positively charged clusters. While the cluster radius is — apart from a 
small spill-out — well-defined for big clusters this is not the case for small clusters. In 
the following we will use (jlj) down to the single, isolated atom (where R = rg) and hence 
plot the respective binding energy vs l/rg. Even with the uncertainty in the definition 
of the single atom or ion radius kept in mind, the binding energy for the single atom or 
ion in general does not lie on the metal sphere line predicted by Hence, there is 

a transition region (hatched area in Fig. |2]) in which the binding energy is expected to 
deviate from the metal sphere line towards the single atom or ion result (filled colored 
circles in Fig. |2]). 

§ Calculating the work needed to move the electron from the metal sphere surface to infinity using the 
image charge-method yields a formula different from ([T4| . Only if the work is taken into account that 
is required to move the electron from inside the metal sphere through the surface, (fT4| is recovered. 
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3. Results and discussion 

Table [H summarizes the PZ-SIC xLSD results for the smallest and largest system we 
consider in this work, namely the single, isolated Na atom {A = 1), and the Nasg 
cluster {A = 58). One sees that in the single atom-case the removal of the 3st valence 
electron and subsequent relaxation gives Ai^ss-i^^rixd = 0.191, in good agreement with the 
experimental ionization potential 0.18!ijj]. The same procedure for a 2pJ, yields a binding 
energy Ai?2p4,rixd = 1.395, i.e., ~ 38 eV, in very good agreement with the experimental 
result for the '^P2 XPS peak [23^. The removal of a 2p'|"-electron gives the slightly 
different and higher binding energy A£'2p-|~,rixd = 1.403. Looking at the KS orbital 
energies one recognizes that — ess^^rixd = 0.186 is slightly below A^Ess^rixd = 0.191. 
Hence, Koopmans' theorem is almost fulfilled for the valence electron, which is not 
surprising because relaxation effects are expected to be small if a loosely bound outer 
electron is removed. 

The difference between KS orbital energy and relaxed energy difference is larger for 
the 2p-electron removal: —e2pi = 1.452 > A_E'2p4.,rixd = 1.395. Obviously, the relaxation 
energy does not fully cancel the "non-Koopmans energy" [llj. 

For Na58 we see in the right Table [U that the binding energy of the valence electron 
^-^'ig.rixd is reduced compared to the atomic expected from the metal sphere 

model for an initially neutral cluster. The binding energy Ai?2p,rixd = 1.308 also shifts 
to lower values with increasing cluster size. However, the absolute value of the KS 
orbital energy increased with respect to the atomic value. Hence, naively estimating 
the photoelectron peak position in an XPS spectrum to be ekin = e2p + ^ is not only 
inaccurate but even may predict peak shifts into the wrong direction. 

3.1. Energies vs 1/R 

Fi gUre[3j3 shows Ail^valence.rlxd) '^-£'valence,unrlxd5 and |evalence| VS 1/R. Wc obsCrVC the WcU- 

known zig-zag behavior of the binding energy because of shell-closures [13], EU [25] before 
the metal-sphere result is approached for increasing cluster size. For the removal of a 
valence electron the KS orbital energies |evaience| are very close to A-Evaience.rixd, and both 
develop towards the metal sphere behavior as l/i? decreases. The relaxation energy 

^valence -^valence, unrlxd ^-^valence.rlxd is Very Small. 

The typical valence electron result for binding energies in clusters shown in Fig. [3b 
is known for decades [T3[ [2H [25] . The qualitatively new aspect coming into play when 
core-electrons are removed is the reorientation of the AZ + — 1 electrons after the 
emission of the core-shell photoelectron. There are (at least) three options: 

(i) the photoelectron is slow enough to gain the relaxation energy which is set free due 
to the screening of the core-hole by the other electrons but fast enough to escape 
before the core-hole itself decays; 

II A pure xLSD calculation yields only 0.179, showing that the PZ-SIC greatly improves the result. 
% Because spin-orbit coupling is neglected, we do not differentiate fine-structure in our KS calculations. 
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Table 1. PZ-SIC xLSD results for the single, isolated Na atom (left) and the Nasg 
cluster (right). The subscripts a and c on the right hand side indicate whether the level 
belongs to the central atom or jellium cluster, respectively. The Nass cluster is a spin- 
neutral system with closed shells (i.e., A = 58 is a "magic number") so that binding 
energies are independent on whether a spin-up or spin-down electron is removed. 



Nai orbital 


energy 


Nass orbital 


energy 


1st 


-40.333 


ISa 


-40.367 


2st 


-2.425 


2s, 


-2.465 


2pt 


-1.453 


2pa 


-1.493 


3st 


-0.186 


3s^ / Is^ 


-0.240 


Isi 


-40.333 


IPc 


-0.194 


2si 


-2.424 


Ide 


-0.166 


2p; 


-1.452 




-0.171 


Eo 


-162.217 


2pc 


-0.139 




-160.822 


Ifc 


-0.138 


-Esst.rlxd 


-162.027 


Igc 


-0.108 


^-£'2p4,,rlxd 


1.395 


Eo 


-351.215 


A-E'3st,rlxd 


0.191 


-E'2p,rlxd 


-349.906 






-E'lg,rlxd 


-351.101 






^-E'2p,rlxd 


1.308 






AE'ig^rlxd 


0.114 



(ii) the photoelectron escapes so rapidly that it cannot pick up the relaxation energy; 

(iii) the photoelectron is slow enough to gain energy not only from relaxation but also 
from Auger-like core-hole decay. 

Signatures in XPS spectra that would support option (iii) were not observed in 
experiments [T6| ITf] . Instead, slow "secondary" electrons are measured, which are likely 
generated via Auger decay of the core-hole after the photoelectron is already gone. 
Indeed the estimated photoelectron escape times are sub-femtosecond for the cluster 
sizes and photon energies considered in this work and thus faster than Auger decay . 

Option (i) is supported by the experiment because the experimental binding 
energies ([6]) for sufficiently big clusters follow again the metal sphere prediction with 
only the work function shifted to the respective bulk binding energy. 

Option (ii) would manifest itself as an increased binding energy of the photoelectron 
as compared to the situation where all the relaxation energy goes to the photoelectron. 
Instead, the experimentally measured binding energies in [161 HZ] drop below the metal- 
sphere result, which ultimately must be so because the single ion value (Pb~ in [161 IE]) 
lies below the metal-sphere line. Hence, one has to distinguish (at least) two effects 
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Figure 3. Relaxed energy-differences (black □), unrelaxed energy-differences (blue 
A), and absolute value of KS orbital energies (green 0) vs inverse cluster radius 
= {A-^/^Ts)^-^ for the removal of an atomic 2p electron (a) and a cluster valence 
electron (b). The small annotations close to the symbols in (b) indicate the number 
of atoms A in the cluster. Red crosses indicate experimental data points, red lines 
[several dashed in (a) and one solid in (b)] the metal sphere-model prediction ([T4l) for 
Ze = 0- For the valence case (b) there is almost no difference between relaxed binding 
energy, unrelaxed binding energy, and the KS orbital energy. However, for the removal 
of a 2p core-shell electron the differences are large and increase with increasing cluster 
size. The difference between unrelaxed and relaxed binding energy is the relaxation 
energy, as indicated by the curved arrow in (a) . The purple -I- symbols in (a) show the 
minimum absolute value of the 2p KS orbital energy for the relaxed configuration as 
the 2p ground state fractional occupation number is varied between 1 and (see text 
for discussion). 



here: the abihty of the other electrons to screen and the abihty of the photoelectron to 
pick up the relaxation energy due to screening. The first one seems to be relevant in 
the experiments fl6\ [T7] so that indeed the metal-to-nonmetal transition in Pb^ around 
A = 20 could be probed. The situation for initially neutral Na clusters is quantitatively 
different. In fact, we find — apart from the shell-oscillations — little deviation from the 
metal-sphere behavior down to the single-atom limit in Fig. |3}d, in agreement with 
experiment [27] . 

Figure [3^ shows the PZ-SIC xLSD result for the removal of a 2p core-shell electron 
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from the central, all-electron Na atom. Again, the relaxed energy difference follows 
approximately the trend predicted by the metal sphere model, although with substantial 
fluctuations on the eV-level (leV ~ 0.037 a.u.). These fluctuations are also correlated 
with the filling of electronic shells but in a less obvious way as compared to the valence 
case where maxima in the (absolute value of the) binding energy occur at half-filled 
and full shells. It would be desirable to move towards bigger cluster in order to see 
whether our model really approaches the metal-sphere result for the screened 2p core- 
hole. However, it is numerically very demanding to proceed significantly towards the 
bulk limit. Note that 1/R = 0.05 already amounts to A = 125. 

Figure clearly shows the importance of the relaxation energy = A£'2p,unrixd — 
A-£'2p,rixd (i-e., the difference between blue triangles and black squares). As expected, 
the relaxation energy increases with the cluster size. For the biggest cluster treated in 
this work, Nass, we have T,^p^ = 0.3, i.e., 8.2 eV. 

The KS 2p-orbital energies |e2p| of the ground state (green diamonds) are between 
the relaxed and the unrelaxed energies, as it is expected from PZ-SIC xLSD [H]. We 
checked how the KS orbital energy shifts as the fractional occupation number /2p is 
decreased from 1 to in the always re/axed system. We observed that |e2p(/2p)| decreases 
below A_E'2p,rixd (purple + symbols in Fig. [3^) and then, for 0.2 < /2p < proceeds 
towards the value for the core-shell-ionized system. If the unknown exact xc potential 
could be employed the orbital energy would jump discontinuously as /2p changes from 
< £ < 1 to 0. This is the so-called "derivative discontinuity" [28j. The PZ-SIC 
xLSD mimics this behavior but smoothes-out the discontinuity. We see that the metal 
sphere result in Fig. [3^ lies between the relaxed binding energy and the minimum KS 
2p-orbital energies min |e2p(/2p)|- Fractional occupation numbers are related to Slater's 
transition state theory where the idea is to consider a (relaxed) configuration "half way" 
between initial and final state. Indeed, e2p(/2p — 0.5) + JtjjJ would predict a photoelectron 
peak in XPS spectra close to what is expected from the metal-sphere model whereas 
e2p(/2p = I) + hu would be completely off and — apart from fluctuations — erroneously 
increasing with increasing cluster size. 

Concerning the metal-sphere result for the 2p core-level binding energy, it should 
be noted that the 2p bulk limit is not precisely known. Moreover, the experimental 
peaks are typically splitted over ~ 0.4 eV because of flne-structure. We estimated the 
2p bulk value from extrapolating the binding energy obtained from the large-cluster 
XPS spectrum in [23]. As a consequence, the core- level metal sphere result in Fig. [3^ 
may have an uncertainty of ±1 eV ~ ±0.04 atomic units (red dashed area). 

The experimental results for the core-level binding energies as a function of the 
cluster size [16l [T7j follow the metal-sphere prediction down to a certain minimum 
cluster size and then deviate, ultimately approaching the single atom or ion result, 
which in the experiment with Pb~ clusters was below the metal sphere binding energy 
for the 5d and 4f core-shells (unlike for the 2p shell in Na above). A too rapid escape 
of the photoelectron would lead to an increased binding energy. It thus seems that 
the relaxation energy that is set free due to screening goes indeed to the photoelectron 
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and that the deviation from the metal sphere-model is an element-specific electronic 
structure-effect. In fact, it is expected that from a certain minimal cluster size on the 
spherical jellium approximation breaks down. However, our extended jellium model has 
the exact single atom-limit, so it may deviate at a too small cluster size from the metal 
sphere line but it ultimately will. The only exception arises when the single atom or 
ion binding energy itself lies — perhaps by chance — on the metal sphere line. Indeed, for 
sodium this is almost the case, as is seen in Fig. [31 

3.2. TDDFT study of the screening dynamics 

In order to understand why the relaxation energy that is set free due to screening goes 
to the photoelectron we determined the relevant time scales. To that end we removed at 
time t = instantaneously a 2pJ, electron from the central, all-electron Na atom in the 
Nasg-cluster, and propagated the KS orbitals thereafter. Figure H] shows the difference 
in the radial spin densities 



as a function of time t and radial position r. The spin-down density shows a very 
similar behavior and is therefore not shown. At small r < 1 mainly the charge density of 
the inner electrons (up to the 2p-shell from which the photoelectron has been removed) 
screens on a very fast time-scale (~ 2 a.u. ~ 48 as) and continues to oscillate. This 
is simple to understand: after the sudden creation of the core-hole, the KS orbitals 
experience a new potential in which they are not eigenstates anymore. If the potential 
was stationary one would expand the KS orbitals in the new eigenstates, and the 
oscillation periods would be given by the inverse eigenenergy-separations. In TDDFT 
the situation is more involved because the potential is nonlinear, i.e., it depends on the 
KS orbital densities. Moreover, we cannot employ linear response theory here because 
the change in the fractional occupation f2pi number is not small but unity. Nevertheless, 
the energy width spanned by the inner-electronic levels is such that it yields the screening 
time-scale visible in Fig. HI As can be inferred from the density oscillations for bigger 
r, the valence electrons screen on a time-scale of ~ 20 a.u. ^ 0.5 fs, i.e. still sub-fs, 
in agreement with Ref . [29] Q The decay of the density oscillations seen in Fig. H] is 
probably underestimated because spontaneous emission, electron ion collisions, energy 
transfer to ionic degrees of freedom, and perhaps other possible dissipation channels are 
absent in our TDDFT approach. 

The time- dependent calculations clearly reveal that the screening dynamics, and 
thus the relaxation, are fast enough to provide the escaping photoelectron with the 
relaxation energy. For instance, the kinetic energy of a 2p electron in a Nass-cluster 
that is removed by a 40 eV photon is ~ 4.6 eV, and it takes ~ 0.6 fs to escape from 
the cluster. Hence, it certainly has enough time to pick up the relaxation energy from 



+ In Ref. [29]. the dynamics after the sudden immersion of a negative charge in a jelhum cluster was 
investigated. 




(17) 
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Figure 4. (a) Contour plot of the difference in the radial spin density 47rr^[ri''~(r, — 
"Ipj. rixd('')] N^^-ss vs time t and radius r (both logarithmically) after the sudden 
removal of a 2p4, electron at time t — Q. (b) Line-outs for r = 0.51 (black), r — 2.51 
(orange), r — 14.01 (blue), showing the different time-scales involved, i.e., "inner- 
atomic" and "inter-atomic" , respectively, as indicated in (a) . 

its own shell (and the shells below), and even the time-scale for the valence electron 
screening is comparable. 

4. Conclusions 

A jellium model with a central all-electron atom has been introduced in order to model 
recent experiments on core-shell ionization of metal clusters as a function of the cluster 
size. Density functional theory in exchange-only local spin-density approximation with 
a Perdew-Zunger self-interaction correction has been employed to calculate the self- 
consistent electron configurations of the model applied to sodium clusters. The relaxed 
binding energy follows the trend predicted by the metal sphere model and confirmed 
experimentally. However, shell-structure-related fluctuations more pronounced than 
those for valence electron-removal are observed. Both the unrelaxed binding energy and 
the initial, ground state Kohn-Sham orbital energy are completely off, demonstrating a 
substantial relaxation energy, the importance of screening, and the irrelevance of initial- 
state core-shell Kohn-Sham eigenvalues. Predicting a photoelectron peak by simply 
adding the photon energy to the initial Kohn-Sham core-level energy is, in general, 
doomed to fail. Instead, for the exchange-correlation potential used in this work Kohn- 
Sham eigenvalues evaluated for fractional occupation numbers (in the spirit of Slater's 
transition state theory) are closer to the relaxed binding energy and metal-sphere result. 
Time- dependent density functional theory has been applied in order to investigate 
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the screening dynamics after the sudden removal of a Na 2p core-electron. Inner-atomic 
screening due to electrons in the same and other core-shells of Na occurs within a few tens 
of attoseconds. Inter-atomic screening by the cluster valence electrons occurs rather on 
the hundreds-of-attoseconds time scale. The inner-atomic time scale is, of course, core- 
shell dependent and may be even faster for deeper-lying shells and heavier elements while 
the screening time scale in metal clusters due to valence electrons is less sensitive and 
depends mostly on the number of valence electrons per atom. Whether the photoelectron 
"has enough time" to pick up the relaxation energy depends on the ratio of escape time 
to relaxation time. The higher the photon energy and the smaller the cluster the higher 
should be the net binding energy. However, the deviation from the metal-sphere model 
observed so far experimentally is most likely dominated by electronic-structure effects. 
In order to verify this assertion future studies will go beyond the jellium model. 
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